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Abstract. We study the dislocation formation in 2D nanoscopic islands with two methods, 
, the Molecular Static method and the Phase Field Crystal method. It is found that both methods 

. indicate the same qualitative stages of the nucleation process. The dislocations nucleate at the 

Q ' film-substrate contact point and the energy decreases monotonously when the dislocations are 

O . farther away from the island-wetting film contact points than the distance of the highest energy 

barrier. 



m 

' 1. Introduction 

The shape and size of islands resulting from growth processes has been a subject of numerous 



o 



CN ; recent experimental and theoretical studies Pl^lSllElElElElEllIIllTSllTllllSj. In particular, 

■ when the growing islands are of nanoscopic size the central issue is the possible spontaneous 
i self-organization of islands into arrays of islands with a narrow size distribution. Such cases 

offer immediate technological applications in modern nanotechnology. The details of such self- 
. organization processes are poorly understood, however. In particular, there are still uncertainties 

r> I as to whether the observed shapes and sizes of growing islands in heteroepitaxy correspond to 

■ thermodynamic equilibrium state of minimum free energy or they are limited by kinetic effects. 
In the earlier analytical and numerical studies, in addition to assuming predefined shapes for 
the islands, the role of dislocation nucleation has not been included. 

In the present work we study the dislocation nucleation in two dimensional islands with two 
methods, the Molecular Static method and the Phase Field Crystal method (PFC) ITOj. The 
appeal of the latter method is that we do not need to initialize the transition path by hand. 
The way in which the two models reach the timescale of the nucleation process is quite different 
and our focus is on the qualitative stages of the process and on the agreement between the two 
approaches. 

In the first section on this paper we explain the methods. We use the first, atomistic, model to 
study the strain relief caused by dislocation formation by minimizing a transition path between 




Figure 1. The minimum energy path be- 
tween coherent and dislocated configura- 
tions is obtained with Nudged Elastic Band 

method in system with a 411- atom island 

"reaction coordinate and 5 layers thick wetting film [l9]. The 

dislocations are nucleated at the island- 

.S i av!y!> wetting film contact points. At a certain 

■vaS^^ distance from the contact point the energy 

:!?!?!?!?^^^ ^ ^^^^ barrier is highest. After this point the en- 

:::::s^^ ergy is lowered monotonously as the dislo- 

^rrt - v— r'^ cations move to the island centre. 



the coherent and dislocated states with Nudged Elastic Band (NEB) Method [20]. Then we 
explore the same process with the PFC method in two dimensions. The nucleation is seen to be 
qualitatively similar with both approaches. 

2. Models and Methods 

2.1. Molecular Static Method 

We adopt the 2D model used previously in Ref. [19], Interactions between all atoms in the system 
are described by a modified Lennard-Jones (LJ) pair potential [16] V{r) with two parameters, 
namely the dissociation energy ^ab and the atomic equilibrium distance Tab, 
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with a smooth cut-off as in [16j. The indices ab of the equilibrium distance are ss, ff and fs for the 
substrate-substrate, adsorbate-adsorbate and adsorbate-substrate interactions, respectively. We 
model the attraction and lattice mismatch between the two materials by changing the potential 
depth ffab and the equilibrium interatomic distance Tab between the substrate and the film. 
The substrate-substrate interaction parameters were set to ess = 3410.1 K and r^^ = 2.5478 
A corresponding to a Cu substrate [16]. In terms of the lattice mismatch rg = (1 + /)rss and 
Tfs = (rgf + rss)/2. A positive mismatch / > corresponds to compressive strain and negative 
/ < to tensile strain when the adsorbate island is coherent with the substrate. Analogously we 
define an interaction parameter k = {£ss—£sf)/£ss which describes the relative difference between 
the potential minimum depths in the substrate-adsorbate and absorbate-adsorbate interactions. 
A negative value of k < corresponds to an effectively attractive and positive «; > to a 
repulsive substrate. The potential depths can be written as eg = Ess, £fs = (1 + ^)£ss for 
the adsorbate-adsorbate and substrate-adsorbate interactions, respectively. We used periodic 
boundary conditions in the horizontal direction. Two bottom layers of the substrate were held 
fixed to simulate a semi-infinite substrate while all other layers were free to relax. 

2.2. Phase Field Crystal Method 

The phase field crystal method (PFC) has been applied in a variety of problems related to 
non-ideal crystal structures [TU]. The method finds the ground state of a given truncated free 
energy functional without any reference to the instantaneous position of individual atoms. This 
enables us to find the dislocated ground states even if they were behind kinetic barriers, which 
would be insurmountable for atomistic simulations. The free energy is 
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Figure 2. The figures in the upper row are PFC 
configurations at different stages of dislocation 
, formation. The figure in the left shows the free 

*® energy landscape as a function of the image 

number in arbitrary units. The figures in the 
upper row correspond to images 1, 5 and 8. The 
misfit / = —7.8% and the pinning potential 
\ depth is half of the amplitude in the periodic 

_^ , , ^ phase. The line in the energy landscape is drawn 

2 4 6 8 to guide the eye. The dislocations positions are 

IMAGE indicated by the square loops. 



and the parameters '(9 = 1/2 and q = A;(l+//100), while k is the period of the pinning potential to 
be defined later. The components of the vector r = (x, z). We use periodic boundary conditions 
in the horizontal and mirror boundary conditions in the vertical direction. The function /i is a 
pinning potential, which zero, if 2z > 3\/^TT/k and otherwise 

h{r) = h cos{kx)cos{kz/V^) -cos{2kz/Vs)/2 , (3) 

where k = b^/3/2. The factor 6 ~ 1 ensures that the potential is continuous in the horizontal 
direction. The amplitude h is the pinning potential depth and it plays the same role as the 
interaction parameter in the Molecular Static method. The difference of the pinning potential 
and triangular phase periods k and q is adjusted to give the wanted lattice mismatch. The time 
step At = 0.005 and the grid size Ax = 0.7854. The simulation was done on 1024Ax x 512Ax 
grid. The average density g = 0.37. This corresponds to a point in the middle of the constant 
phase-triangular phase coexistence in the one mode approximation of the phase diagram of Eq. 
[21 The initial state has an island of half-disc shaped region of triangular phase on the bottom 
edge of the simulation cell, see Figj2j This island is surrounded by constant phase. The initial 
state is evolved according to 

where C is Gaussian distributed conserved noise [lOj . 



3. Results and Discussion 

Fig. [1] shows an energy path between coherent and dislocated configurations [19]. The misfit 
/ = 8% and k = 0. The initial chain of states is obtained by simple linear interpolation. This 
chain is minimized with NEB method |20]. The results show an energy barrier for the initial 
detachment of the dislocation from the island-wetting film contact point. When the distance 
to the contact is larger than the distance of the highest energy peak, the energy decreases 
monotonously as the dislocations approach the island center. Finally the defects stop at a finite 
distance from each other. 



We used the PFC model to make a similar calculation. We used mismatches of ±7.8% and 
pinning potential depth 0.5 times the amplitude of the triangular phase. The initial state is the 
upper half of a disk standing on the simulation cell edge where the pinning potential is effective. 
For compressive mismatch, the island does not form dislocations but adopts a concave shape 
with reduced size of the island-wetting film contact area. 

In the tensile case dislocations are nucleated symmetrically from both sides of the island 
base. This is shown in Fig. [21 The defects move towards the island centre. The free energy 
landscape is given in arbitrary units as a function of the image number. It shows the combined 
effect of dislocation nucleation and facet formation. These effects could be studied separately 
by starting from dislocated or faceted initial configuration. The general shape of the curve in 
Fig. [2] is similar to that of Fig. [TJ 

The island shape is almost rounded in the compressive and strongly faceted in the tensile case. 
The asymmetry of tensile and compressive has been studied also in [17\ I18j. This asymmetry is 
clearly visible in the PFC model, because it does not nucleate dislocations in the latter case. 

4. Conclusions 

We have studied dislocation nucleation in 2D strained heteroepitaxial nanoislands. We have 
employed both an atomistic model and phase field crystal method. The results are qualitatively 
consistent with each other. The dislocations nucleate near the island-wetting film contact under 
tensile strain. After an initial energy barrier the dislocation detaches from the contact point 
and moves to the middle of the island. During the move the energy is lowered gradually. The 
island shapes in the tensile case were observed to be more faceted than in the compressive case. 
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